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The semi-inclusive cross section for producing slow protons in charged current deep inelastic 
^ ■ (anti-) neutrino scattering on protons and neutrons is calculated as a function of the Bjorken x 
^^ • and the proton momentum. The standard hadronization models based upon the colour neutral- 
-(— > 



ization mechanism appear to underestimate the rate of slow proton production on hydrogen. 

The presence of the virtual mesons (pions) in the nucleon leads to an additional mechanism 
^ I for proton production, referred to as spectator process. It is found that at low proton mo- 
^ ■ menta both mechanisms compete, whereas the spectator mechanism dominates at very small 

momenta, while the color neutralization mechanism dominates at momenta larger than 1-2 
^ . GeV/c. The results of the calculations are compared with the CERN bubble chamber (BEBC) 
H i data. The spectator model predicts a sharp increase of the semi-inclusive cross section at small 

X due to the sea quarks in virtual mesons. 



PACS numbers: 25.30.Pt, 13.85.Hi 



1 Introduction 

Inclusive deep inelastic lepton scattering off nucleons is a well established tool for investigating 
the parton distributions ||l|. Hadrons, which are not measured in inclusive experiments, are 
produced mainly in the colour neutralization process; the string models being the state of art. 
It is expected that the measurement of final state hadrons will lead to a better understanding 
of the mechanism of hadron production as it contains information about hadronization process. 

In semi-inclusive deep inelastic reaction (SIDIS) one observes one or more hadrons in co- 
incidence with the scattered lepton. In the simplest approximation the semi-inclusive cross 
section is assumed to factorize into an inclusive cross section (proportional to the deep-inelastic 
structure functions) and a fragmentation function |0, |[. 

In the past the study of semi-inclusive deep-inelastic lepton scattering has been restricted 
mainly to the detection of high energy hadrons, which originate from the fragmentation of the 
leading (struck) quark. The production of slow protons in coincidence with muons has been 
studied in the (anti)neutrino induced charged current reactions at CERN @, |^. In a more 
recent analysis of the CERN experiment WA21 [^ the fraction of slow {p < 0.6 GeV/c) protons 
was determined. In Ref. the rate from Ref. [^] was reproduced with a constant fragmentation 
function, assuming that only protons (no other baryons) are produced. The rate of slow proton 
production strongly depends on the behaviour of the fragmentation function at the kinematical 
limit, i.e. the light cone momentum fraction z = p~ /m close to 1. The counting rules suggest 
for the diquark fragmentation function rather a linear, D2q{z) = A{1 — z), than a constant 
dependence on z. Furthermore, other baryons (neutrons, isobars, etc.) can be produced which 
leads to a suppression of the relative proton yield. 

In order to illustrate this point we show in Table 1 the relative yield of protons and neu- 
trons including both direct and sequential (intermediate A) production obtained from a more 
realistic, yet simple, quark-diquark model ^. For comparison we give the number of pro- 
tons/neutrons per event obtained from the Lund string model 0. Both arguments discussed 
above suggest that the agreement with experimental rate in Ref. P| has to be reconsidered. It 
will be discussed later that a more realistic calculation leads to a deficit of slow protons com- 
pared to the experimental rate [Q . The violation of the Gottfried sum rule observed by NMC 
|T0| strongly suggests the flavour asymmetry of the sea of light quarks in the proton. This can 
be naturally accounted for by the presence of the pion (meson) cloud in the nucleon |ll|] - |jT8| 



The pion cloud model has also been used successfully in inclusive DIS to explain the antiquark 
distributions [l^, |l^ . 



The presence of the meson cloud may lead to an additional mechanism for the slow proton 
production in deep inelastic lepton scattering pO, pTl 521 Bl]. For instance, in the (anti)neutrino 



induced reaction the virtual boson W smashes the virtual pion (meson) into debris and the 
nucleon (baryon) is produced as a spectator of the reaction (see Fig.l). The presence of the 
IT A Fock component in the nucleon leads to the production of a spectator A which decays into 
a pion and nucleon. 

Up to now no fully consistent treatment of both the standard fragmentation mechanism 
and the meson exchange mechanism exists. We note that it is inconsistent to add the cross 
sections for the two processes. In the present paper we discuss in a simple hybrid model how 
the two mechanisms can be taken into account in a consistent manner which ensures baryon 



number conservation. It is also not fully consistent to replace the four-quark fragmentation by 
the meson exchange mechanism due to the following two reasons. First of all, as discussed in 
Ref. [0, the virtual pions (mesons) give rise both to the sea and valence quark distributions. 
Secondly ||18|, ^ the sea generated by mesons constitute only a fraction of the total sea. 



It is the aim of this paper to quantitatively predict the contribution of the meson exchange 
mechanism to the slow proton production rate in (anti-)neutrino charged current deep inelastic 
scattering from nucleons. In addition the momentum distribution of protons is calculated and 
compared with existing experimental data 0]. Finally the x-dependence of the coincidence cross 
section for the experimental momentum bins is calculated and compared with experimental data 
fromRef.e. 



2 Deep inelastic (anti-) neutrino scattering 

2.1 Inclusive scattering 

In leading order electroweak theory the inclusive, differential, charged current (cc) cross section 
for the neutrino/antineutrino scattering {h{i,i')X) can be expressed in terms of the structure 
functions -Fi, F2, F3 

'^''''^ - ^^ ■[y'xF,ix,Q') + il-y)F,ix,Q')±iy-y^)xF,ix,Q% (1) 



dxdydQ^ ^(1 + ^)^" .V ,- . ■ V .. .V ,^ . V. 2 
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Here Gp = 7r/(\/2 sin^ 6'h/M^) is the Fermi coupling constant, Mw is the 1^-boson mass and 
y = {Ep — El) /El. Assuming the Callan-Gross relation, 2xFi = F2 and the standard parton 
model the y integrated cross section is given by 

^ ^^^ - !l£ xT^ (x O'^) (2) 

dxdQ^ 7r(l + QyM^f^^'^''' "^ ^ ' ^^^ 

where J-'^d^, Q"^) can be expressed in terms of the quark distributions 

^:,(x,Q2) = d{x,Q^)+u{x,Q^) + ]^{-y+y-)[d{x,Q^) + M{x,Q% 

Tl{x,Q^) = d{x,Q^)+u{x,Q^) + ^{-y + y^)[d{x,Q^) + 3u{x,Q^)] (3) 

for neutrino and antineutrino induced reactions, respectively. In the expressions above the 
distributions of heavier quarks have been neglected. In the scaling limit and transferred four- 
momenta Q^ <^ M^ the Q^-integrated cross section is 

^ = — ^ — ^-^^^' ^^^ 

where 

g^^^x) =2x[d{x)+u{x)/3], 

Glix) = 2x[d{x) + u{x)/3]. (5) 

This formula will be useful in the further analysis. 

2.2 Semi-inclusive scattering: fragmentation approach 

Semi-inclusive charged current (anti) neutrino reactions, where a slow proton is measured in 
coincidence with the muon, are of special interest. The Q^-integrated cross section for the 
production of a baryon B (nucleon, delta) can be written as 

da[N{i,i'B)X] _ G^rriME^,^^ 
dxdzdp\ 71 

Uvr''\-)DS'i~^yT) + Sf^^''\x)D!Jfizy,)]/il - X). (6) 
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Here Vf ' (x) and Sr (x) are functions proportional to the valence and sea quark distri- 



'/ 



butions, respectively, with coefficients being reaction and flavour (/) dependent, D^e^T-isi^ , Pt) 
is the target debris -^ proton fragmentation function, with z = p~ /rri]^, z = z/{\ — x), and px 
is the transverse momentum of the emitted proton (baryon). Usually a factorized form of the 
fragmentation function is assumed 



Dl,^,,r^,pl)=D{i)^{pl). 



(7) 



Then, neglecting a small production of the BB pairs, the baryon number conservation requires 
for each flavour the following sum rules (assuming only one baryon is produced in the target 
fragmentation region) 
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z) dz = 1, 



(8) 



separately for diquark, D2q{z), and four-quark, D^g^z), fragmentation functions. 

In more advanced string models [^, ^, |2^, ^, ^ the semi-inclusive cross sections can be 
calculated with the help of Monte Carlo methods PO, ^^, 521 H . 



2.3 Pionic contribution to the semi-inclusive cross section 

In the simplest models the nucleon is treated as a system of three quarks. The pion-nucleon 
interaction leads naturally to an admixture of a iiN Fock component in the physical nucleon. 
In the simplest approximation, the Fock state decomposition of the hybrid proton (neutron) 
reads 



\P) 



cos6'|(po)) + sin6' 



sm(f)(,/2/3\n+n) - ^/l/3\n^p) 



-f-cos 



4> (^yi72|7r-A++) - yi73|7rOA+) + yi/elTr+A") 



(9) 



sin0 [yi/SlTT^n) - ^2/3^"^) 



\n) = cos 6^1 (no)) + sin 6^ 

+ COS0 fv^k^^") - \/V3|vr^'A") + Jl/6|7r-A+) 



(10) 



where only nN and vrA Fock states are taken into account. 

The two mixing angles are related to the number of pions (n^ = n^Ar + n^rA) in the pionic 
cloud 

sin^6' = n7r/(l + n^) 
sin^ (p = n.„M /ut, . 

The Fock expansion in Eqs. (^ |10[) leads to the processes which can be classified as spectator 
mechanism (see Fig.l), core fragmentation, and intermediate nucleon (and intermediate A) 
fragmentation with the pion produced as a spectator, respectively. 



Schematically the total cross section can be expressed as 

A 

+ Z[a^{N -^N') + J2^'^iN ^ A^N')] 

A 

N" N",A 

+ P.aE /ac^'(A ^N')+Y: fAcr\A ^ A' ^ N)] . (11) 
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In the formula above ct^p and a^ are cross sections for the spectator and colour neutralization 
mechanisms, respectively. The factors Z = cos^ 9, P^n, P-kA are probabilities of the core, and 
irN and vrA Fock components. The f^r and /a factors are the relative weights of the nucleons 
(deltas) in the meson cloud model as dictated by the isospin symmetry. A fully consistent cal- 
culation would require to include the momentum distribution of virtual baryons in the nucleon, 
which goes beyond the scope of the present paper. In the actual calculation we have taken 
into account only the first two lines in Eg. (|Tl|) . We note that the nucleon can be produced 
in the direct way (M = 7r,B = N) or as a two step process: direct production of isobar A 
(M = 7r,B = A) and subsequent decay of the resonance {A ^ N + it). 

In calculating a^ (see Eq.(|ll])) one should bear in mind that part of the the sea quark 
distributions in bare nucleons (deltas) is already taken into account by the pionic component. 
Therefore in Eq.(^) the sea structure functions Sf{x) should be replaced by the difference 
between the experimental ones and the ones computed in the pion spectator approach. In 
practice we find that at Q^ = 4 (GeV/c)^ this is equivalent with a reduction of the four quark 
fragmentation by a factor 0.5 — 0.7 |l^, Q. 

The probability density to find a meson with longitudinal momentum fraction t/M and four- 
momentum of the virtual meson squared, t, (or alternatively transverse momentum p^), referred 
to as splitting function, quantifies the presence of virtual mesons in the nucleon. The splitting 
function f{yM,t) for the first process is 

fMVM.t) - ^^yu ^^_^,^, . (12) 

The splitting function for the ttA Fock state is 

where M+ = m/\ + uin and M_ = mA — thn ■ The F{yM,t) are vertex form factors, which 
account for the extended nature of hadrons involved. The form factors used in meson exchange 



models are usually taken to be a function of t only. As discussed in Refs. [jT5|, ^ such form 
factors are a source of momentum sum rule violation. 

The pion (meson) emitted by the nucleon can interact with a virtual IV-boson. In general 
the cross section for the semi-inclusive spectator N{i,i'B)X process depicted in Fig.l can be 
written as 

da'P{N{i,i'B)X) ^ . 2^da'^' / / ^ t^^^ 



Integrating over unmeasured quantities one recovers the known one-pion exchange contribution 
to the inclusive cross section 

dyAifniyM) 1 , (15) 



dx Jx dx 

where the ^^ir ^^ ^^^ Q^-integrated cross section for the deep inelastic scattering of i from the 
virtual pion. In practical calculations the on-mass-shell in cross section given by Eq. (|) can be 
used. 

The probability to find a nucleon A^' associated with the lepton £' at given value of the 
Bjorken x can easily by calculated as 

P«y(x) = I < itihflh^ > \'J^ dyMUN'/N{yM)Q'Z{^/yM) 

+ E I < 1^3, \t^\\t^ > \' f dyMUA/N{yM)0'.:{x/yM)P{A ^ ttN'). (16) 

.A A A Jx 

'■3 

The first term in Eq.(|TBp describes the direct nucleon production whereas the second term 
corresponds to the two step mechanism. The Clebsch-Gordan coefficients account for the 
relative yield of a given isospin channel. The sum involves all charge combinations of the pion 
and A which lead to the final nucleon N'. The probability of the A resonance decay into the 
channel of interest (A^') can be calculated assuming isospin symmetry 

PiA^nN') = \<ltllt^'\l4>\'. (17) 



The formula given by Eq. ([Tq ) can be easily generalized to the case with the proton (nucleon) 
in a given momentum range (pmin < p < Pmax)- In this case the following substitutions are 
necessary 



fnN'/NiyAl) -^ I fiTN'/N{t-,yM)dt , 

min 
/■*max 

f7TA/N{yM) -^ J^ f7rA/Nit,yM)Pit,Pmm,Pmi,x)dt . (18) 

mill 

For the direct process the experimental limits on the proton momentum directly correspond to 
limits for t: 



^min = 2[m^-mjvVSar + ^], 



^max = 2[m^-m^^p^i„ + m^], (19) 

with an extra absolute limit t < — m^y|.^/(l — ?/m)- For the proton production through the 
intermediate delta resonance one has to integrate over all possible t and include an extra prob- 
ability function Pit,Pmm,Pmax)- The -P(t,Pmin,Pmax) is the probability that the intermediate 
isobar A produced with momentum 



Plit) 
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ml (20) 



decays into a pion and nucleon A^' in the given momentum range. It can easily be calculated 
assuming a sharp A resonance and its isotropic decay. Technically it is convenient to calculate 
-P(i,Pmm,Pmax) as a difference 

-P(t,Pmm,Pmax) = P(t,PN' < Pmax) " Pit,PN' < Pmin) (21) 

with self-explanatory notation. In this approximation the probability is 

l^[t,PN' < Pmaxj - 7J J- -—-, -717^ > {^^) 

^ PoiEA/rriA -I) ' 

where po, Eq are momentum and energy of the nucleon in the delta rest frame. 

The momentum distribution of the spectator baryon B produced in the iV(£, £'B)X reaction 
which proceeds via scattering of virtual boson (7, W±) on the virtual meson M can generally 
be written as 



dN{p,Q') ^^JdxJ^dyM[l!:fMBMyM,t)dt] ^(^,Q^) 



(23) 



with ti and t2 corresponding to momenta p ± Ap/2, respectively. The function fMB/N{yM,t) 
can be calculated for different partitions of the nucleon into a meson M (pseudoscalar, vector) 
and baryon B (octet, decuplet) p3 |. 



3 Results and Discussion 

Let us first discuss the probability for production of slow protons as obtained from the hadroniza- 
tion models. (With slow protons we mean hereafter unless specified otherwise, protons with 
momenta in the interval (0.15 - 0.60 GeV/c), corresponding to the experimental bins defined 
in Ref. 0). In Table 2 we list the probabilities calculated in (a) the quark-diquark model from 
Ref. [|], (b) the quark-diquark model corrected for the presence of the pionic cloud (to be dis- 
cussed later) and (c) in the Lund string model 0. In the first two cases the results are rather 
sensitive to the form of the fragmentation function D{z). In this calculation the transverse 
momentum distribution function ^{p^) in Eq-(0) is parametrized by a Gaussian form 

$(p^) = ^exp(-p^/a2) (24) 

with 0"^ = 0.3 (GeV/c)^. We find that the use of a constant fragmentation function (as in 
Ref. [^]) yields a result almost consistent with the experimental one |^. The rate obtained 
with quadratic D{z) = 1/2 * 5(1 — z) and especially with the 'triangular' D[z) = 2(1 — z) 
(counting rule motivated for the diquark hadronization) hadronization functions is lower than 
the experimental rate. The quark-diquark model ^ gives only isospin factors for the valence 
quark (diquark) hadronization. 

To show the sensitivity to details of the sea (4-quark) fragmentation two extreme limits 
have been used in which protons being produced either with zero probability (only diquark 
hadronization) or with unit probability (neglecting production of other baryons), respectively. 
As is seen from the results the contribution of the sea hadronization is rather important for the 
production of slow protons. The total rate obtained in both extreme limits is, however, below 
the experimental rate 0. The Lund string model predicts an even lower rate. 

This supports the conjecture that another mechanism plays a role and therefore we now 
turn to a discussion of the role of pion exchange mechanism as sketched in the previous section. 

In all calculations discussed below exponential vertex form factors 

F{yM,PT) = exp[ ^ ] (25) 

have been used in Eqs.(|l2|,|l3]). In Eg . (125|) MMBiyiPr) is the invariant mass of the intermediate 



meson-baryon system. The cut-off parameters used in the present calculation (A^jv = 1.10 GeV 
and AjrA = 0.98 GeV) have been determined from the analysis of the particle spectra for high- 
energy neutron (p(p, n)X) and A ( p{p, A++)X ) production [jl9|. With these cut-off parameters 



the NMC result for the Gottfried sum rule [|T0| which depends sensitively on A, has been 
reproduced [[l^, |l^ . Furthermore the model describes u/d asymmetry extracted recently from 



the Drell-Yan NA51 experiment at CERN p3. Finally the deep-inelastic structure functions 



of pions are taken from Ref . [^ and those for the nucleon from Ref . |^ 



A quantity of interest [§ is the probability for finding a slow proton in a momentum bin 
Ap as a function of x 

P(a;,Ap) = ^^%M(^)-i. (26) 

ax ax 

9 



The quantity calculated in the spectator model is shown in Fig. 2 where in addition we show the 
direct component by the dashed line (the contribution of the intermediate A is a complement 
to the total P{x) given by the solid line) and the contribution from the valence quarks in the 
pion by the dotted line (the sea contribution is again a complement to the solid line). As 
seen from the figure the direct component dominates the cross section in the whole range of 
X. The contribution of the sea in the pion plays an important role only at small Bjorken x. 
Experimentally this region can only be accessible at high beam energies. 

In order to obtain an absolutely normalized "experimental" quantity the following procedure 
has been applied. In Ref . the ratios of normalized ^-distributions with and without secondary 
protons in the momentum range Ap 



dx 



X, Ap)/(r{Ap) 



have been given, where (1 — Ap) refers to the complement of the momentum interval Ap. 
Taken the absolute normalization i?tot = o'{Ap)/a'tot from Table 1 of Ref. f^ at face value one 
can extract the absolutely normalized quantity (p6D as 



The results of the procedure are presented in Fig. 3 for low (0.15 -0.35 GeV/c) and high (0.35 
-0.60 GeV/c) momentum bins. The results of the spectator model (dashed line) and of the 
Lund string model0 (solid line) are shown in Fig. 3. As can be seen from the figure the string 
breaking mechanism underestimates the experimental data. To get more insight into this result 
it is instructive to consider simpler quark-diquark model P] for the flavour dependent diquark 
branching ratio into nucleons and deltas. In this calculation we have taken the QCD motivated 



fragmentation functions ^6\ 



D2g{z) = C~z^l'^{\ - ~z), Dig{z) = const. (29) 

Since no realistic model for the four-quark branching ratio exists, one can only consider lower 
(diquark fragmentation) and upper (maximal proton production in four-quark hadronization) 
estimates, as for the relative yields in Table 2. The results are shown in Fig. 3 as dotted lines. 
The difference between the lower and the upper limit clearly demonstrates a significant sensi- 
tivity to the four-quark hadronization. The presence of the virtual pions leads to a reduction of 



the perturbative sea by about a factor 2 [gj, which would result in the same reduction of the 
four-quark contribution to P{x). Even the upper limit underestimates the experimental data. 
This suggests a presence of an extra competing mechanism. 

For the proton production on the neutron the situation is somewhat more complicated (see 
Fig. 4). Here the result for the Lund string model falls only slightly below the experimental data 
[^ extracted from the (z/, /i~) reaction on the deuteron (Fig.4a,b). The spectator mechanism 
predicts a very similar result, also almost consistent with the experimental data. For the (z7, /z^) 
reaction (Fig.4c,d) the situation is similar, except that here the corresponding cross sections 
are significantly smaller. Similarly as for the proton production on the proton we present 
the lower (diquark fragmentation) and the upper (maximal proton production in four-quark 
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hadronization) limits of the model discuss in section 2.2. Again a large difference between the 
two limits can be observed. The proton production on the neutron , however, differs from that 
on the proton. It is reasonable to expect here rather the lower limit to be more realistic. The 
experimental data seems to support this expectation. 

In order to better understand the range of proton momenta where the spectator mechanism 
plays an important role, it is instructive to calculate the momentum distribution of protons pro- 
duced in the neutrino (anti-neutrino) induced reactions. The proton momentum distributions 
were determined in the early BEBC experiments at CERN [^]. In order to avoid the resonance 
contributions and minimize the effect of the proton's sea, the events have been selected with 
the constraints W > 3 GeV and x > 0.1. 

The cuts are rather important especially for the low-momentum part of the spectrum. In 
the experiment 0] the wide band (anti-)neutrino beam has been used. In principle, the effect 
of the beam spread should be taken into account. It has been checked, however, that the 
results of the Lund string modelQ are not very sensitive to the beam energy. The results of 
the Monte Carlo simulation for three different beam energies (20, 50, 80 GeV) are shown in 
Fig. 5. The total number of protons per event in the momentum range (0.15 < p < 0.60 GeV/c) 
and corresponding to the experimental cuts[Q, calculated in the Lund string model and the 
spectator model is given in Table 3. Also shown are the experimental results extracted from 
Ref.[§]. Similarly as for the experiment from Ref. the string mechanism produces too little 
low energy protons. In this case the spectator mechanism gives almost total missing strength 
(see Table 3). 

In order to get insight into the total proton momentum distribution, it is instructive to 
consider the simpler quark-diquark model [Q for the flavour dependent diquark branching ratio 
into nucleons and deltas. The momentum distribution calculated with the QCD motivated 
fragmentation functions (Eq.(|29|)) is shown in Fig.6. It was assumed here for simplicity that 
the four-quark hadronization leads in all cases to proton production. This will lead to an 
overestimation of this contribution. The diquark and four-quark hadronization contributions 
are shown separately by the dashed and dotted line, respectively. As clearly seen from the 
figure, the four-quark hadronization populates the slow proton region, whereas the diquark 
hadronization is responsible for the production of rather fast protons. 

The momentum distribution of protons produced in the spectator mechanism calculated 
according to Eq. (|23|) is shown in Fig. 7. The solid line includes both protons originating from 
the ir^p Fock component (direct component) and those originating from the decay of spectator 
deltas. The direct component is shown separately by the dashed line. In the neutrino induced 
reaction the sequential mechanism is as important as the direct production of protons, whereas 
in antineutrino induced reaction its relative contribution is much smaller. This effect is a 
consequence of the dependence of underlying coupling constants on the charge of pions and of 
the difference of F2^'''(x) and F2'''{x) structure functions. In order to better understand the 
underlying reaction mechanism we present also in Fig. 7 separately the contribution from the 
valence quarks in the pion (dotted line). The corresponding contribution from the sea is a 
complement to the total momentum distribution. 

In the experiment of ref.[Q not all protons with low momenta ( p < 0.6 GeV/c) were 
identified. In order to make comparison of the model results in the whole momentum range the 
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inclusion of the proton identification efficiency is neccessary. In the present paper the proton 
identification efficiency has been parametrized as 

e(p) = [l + exp[^:^^r" (30) 

with po = 0.86 GeV/c, a = 0.1 GeV/c and n = 0.25. The resuh of the spectator mechanism 
and the result of the quark-di quark model have been multiplied by the efficiency function e{p) 
and are shown in Fig. 8 by the dashed and dotted line, respectively. One clearly sees that the 
colour neutralization mechanism (dotted line) is dominant for the proton production in a broad 
momentum range. The spectator mechanism (dashed line) plays, however, important role at 
very low momenta. Since in our model the pionic cloud constitutes about half of the nucleon 
sea [^] for Q^ of a few GeV^, for consistency the four-quark hadronization contribution (dotted 
line in Fig. 6) must be reduced by a factor of about 0.5. Therefore the sum 

Nip) = e{p) [N^M + 0.5iV4,(p) + N^{p)] (31) 

consistently includes the two mechanisms. The sum (solid line in Fig. 8) almost coincides with 
the original quark-diquark model and describes the experimental data [Q in a resonable way. 

The considered x-distributions P{x) and proton momentum distributions do not give a 
direct evidence for identifying the meson cloud in the nucleon, although they do not contradict 
its existence. Therefore one should look for other possibilities. It has been suggested by some 
authors to consider the production of slow (for the fixed target experiments) pions. Slow pions 
could be created by a mechanism analogous to that considered in the present paper for the 
slow proton production, as a spectator of the reaction (see Fig. 9). The momentum distribution 
of pions can be calculated in a fully analogous way to that of the slow protons (see formula 
( |23| ) ). In this case all kinematical variables associated with the meson M (pion) have to be 
replaced by corresponding variables for the recoil baryon B (nucleon or delta). The momentum 
distributions of pions (vr^, 7r°, vr") produced in the (z^, /i~) and {y^iX^) reactions calculated in 
this way are shown in Fig. 10, separately for Ntt (solid line) and Avr (dashed line) mechanisms. 
For comparison the momentum distribution of pions generated by the Monte Carlo method 0] 
according to the the Lund string model is shown by the solid line with the centered points. As 
seen from the figure the spectator pions constitute only a small fraction of all pions, therefore 
their identification seems to be rather difficult. 

The relative fraction of different pions differs between the Lund string model and the 
spectator model (see Table 4). This is especially visible for the neutrino induced reaction on the 
proton and antineutrino induced reaction on the neutron. The spectator mechanism predicts a 
huge asymmetry between tt"*" and vr^ production. The reason here is very simple. Both vr" in 
the v+p and 7r+ in the V+n reactions cannot be produced in the scattering of the virtual boson 
W^ off the nucleon of the ttN Fock states. On the other hand their production via scattering of 
W^ from the virtual delta of the vrA Fock state is highly suppressed as here the corresponding 
structure functions vanish in the valence quark approximation. The huge difference between 
the predictions of the Lund string model and the spectator mechanism could in principle be 
tested if one limits the analysis only to the pions produced in the backward directions. This 
part of the phase space is expected to be dominated via the spectator mechanism. A possible 
confirmation of the asymmetry in the production of the positively and negatively charged pions 
would be a strong test of the spectator mechanism and the concept of virtual pions in the 
nucleon. 
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4 Conclusions 

In the present paper the production of slow protons produced in (anti)neutrino induced reac- 
tions on proton and neutron has been studied. We have analyzed the total rate, the Bjorken-x 
dependence of the slow proton production as well as the momentum distribution of emitted 
protons. 

The analysis of the rate of slow proton production in the neutrino and anti-neutrino charged 
current reactions shows that "realistic" hadronization models underestimate the cross section 
measured by the bubble chamber Collaboration at CERN [|]. Although the pion exchange 
mechanism leads to a rather small contribution to the total proton production, it plays an 
important role in the production of slow protons and helps us to understand the deficit of 
the slow proton production as predicted by standard hadronization models based on colour 
neutralization mechanism. 

Even these two mechanisms seem to fail to describe the total rate of slow proton production 
from Ref. 1^. Comparison of the total rate as extracted here from earlier CERN measurement 
m, where more severe restrictions on the Bjorken-x and invariant mass of the produced hadrons 
have been imposed, suggests that the 'missing' mechanism is related to the region of small x 
and low invariant mass of the produced hadrons. The considered x-distributions P{x) and 
proton momentum distributions for (anti-)neutrino charged current reactions do not give direct 
evidence for identifying the meson cloud in the nucleon, although they do not contradict its 
existence. Therefore one should look for other possibilities to prove or disprove the whole 
concept of the meson cloud in the nucleon. 

The recent results of deep inelastic electron-proton scattering at HERA, (ref||3^), open such 
a possibility. The observation of a surprisingly large number of events characterized by a large 
rapidity gap, strongly suggests the presence of a (up to now not really identified) proton with 
approximately beam velocity. These so-called diffractive scattering events, which cannot be 
explained by conventional hadronization codes, have up to now primarily been interpreted in 
terms of Pomeron exchange. Since the Pomeron by definition has vacuum quantum numbers 
this approach predicts that only protons are produced in the final state in contrast to the 
mesonic cloud model which leads to a definite prediction for branchings into proton, neutron 
and isobars. Therefore the planned search for neutron spectators will provide more information 
on the mechanism for the production of slow baryons. 
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Tables 



Table 1. Fraction of protons (/p = Np/{Np + Nn)) and neutrons (/„ = Nn/{Np + A^„)) 

produced in the e + p, u + p and V + p reactions. The fractions have been calculated according 

to the quark-diquark model and the Lund string model for the beam energy 50 GeV. 

The total number of protons and neutrons per event in the Lund string model exceeds unity 

which is caused by the explicit production of baryon-antibaryon pairs. 



beam 


p^p 


p ^* n 


n ^- p 


n —>■ n 


quark-diquark model 










/i^ 


0.5473 


0.4527 


0.2160 


0.7840 


V 


0.9259 


0.0741 


0.5000 


0.5000 


V 


0.5000 


0.5000 


0.0741 


0.9259 


Lund string model 










/^^ 


0.60 


0.47 


0.43 


0.64 


V 


0.78 


0.29 


0.56 


0.52 


V 


0.52 


0.54 


0.31 


0.76 
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Table 2. Fraction (in percent) of events with low (0.15 < p < 0.60 GeV/c) momentum 
protons for neutral current charged lepton DIS and charged current neutrino and antineutrino 

reactions. 



final state baryons 


D{z) 


p + fi^ -^ p 


p + u ^ p 


p + u ^ p 


simple quark- diquark 










model 


constant 


3.47-5.88 


6.70-9.75 


2.74-6.57 




quadratic 


2.15-3.89 


4.22-6.42 


1.69-4.45 




triangular 


0.89-1.64 


1.76-2.70 


0.70-1.89 


improved quark- diquark 










model 


constant 


3.43-5.54 


5.72-8.04 


2.75-6.13 




quadratic 


2.13-3.65 


3.59-5.26 


1.70-4.14 




triangular 


0.88-1.54 


1.50-2.21 


0.70-1.75 


Lund string model 










E = 20 GeV 




2.97 


0.89 


1.44 


E = 50 GeV 




4.37 


1.43 


1.95 


E = 80 GeV 




4.60 


1.42 


1.75 


spectator model 










vrA^ 




0.97 


1.23 


1.55 


ttA 




0.29 


0.52 


0.27 


valence 




0.88 


1.33 


1.28 


sea 




0.38 


0.42 


0.53 


sum 




1.25 


1.75 


1.81 


Expil 






10 


8 
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Table 3. Fraction (in percent) of events with protons in the momentum interval 0.15 < p < 

0.60 GeV/c for the neutrino and antineutrino induced reactions. The results from the Lund 

string model has been obtained with extra limitations: x > 0.1 and W > 3 GeV for three 

different beam energies. 



contribution 


p + u -^ p 


p + u ^ p 


spectator mech. 






ttN 


1.23 


1.55 


7lA 


0.52 


0.27 


sum 


1.75 


1.81 


Lund string model 






E = 20 GeV 


0.98 


1.39 


E = 50 GeV 


1.27 


1.79 


E = 80 GeV 


1.17 


1.58 


Exp.i 


4.74 


5.20 



Table 4. Fractions of different pions calculated in the Lund string model 
energy 50 GeV and in the spectator mechanism. 



for the beam 



reaction 


Lund string model 


spectator mechanism 


vr^ 


TT^ 


7r+ 


TT^ 


TT^ 


7r+ 


H^+p 
I' + p 
V + p 


0.29 
0.23 
0.36 


0.36 
0.36 
0.37 


0.35 
0.41 
0.27 


0.25 
0.01 
0.29 


0.38 
0.24 
0.40 


0.37 
0.75 
0.31 


11^ -\-n 
u + n 
V + n 


0.36 
0.28 
0.42 


0.36 
0.36 
0.36 


0.28 
0.35 
0.22 


0.63 
0.29 
0.71 


0.28 
0.41 
0.25 


0.09 
0.31 
0.04 
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Figure Captions 

Fig.l The spectator mechanisin for the production of slow protons in the (anti)neutrino in- 
duced reactions on the nucleon. 

Fig. 2 P{x, Ap) for 0.15 < p < 0.60 GeV/c calculated within the spectator model (solid line). 
The contribution of the direct proton production is shown by the dashed line and the 
contribution from the valence quarks in the pion by the dotted line. 

Fig. 3 P{x, Ap) for p{v,p)X and p{T7,p)X reactions, for low (0.15 < p < 0.35 GeV/c) and high 
(0.35 < p < 0.60 GeV/c) momentum bins. The "experimental" data has been obtained 
by the procedure described in the text. The result from the spectator model is shown by 
the solid and from the Lund string model by the dashed line. In addition the two limits 
discussed in the text are shown by the dotted lines. 

Fig. 4 P(x,Ap) for n{v,p)X and n{v^p)X for low (0.15 < p < 0.35 GeV/c) and high (0.35 
< p < 0.60 GeV/c) momentum bins. The experimental data are taken from Ref.[Q. The 
result from the spectator model is shown by the solid and from the Lund string model by 
the dashed line. In addition the two limits discussed in the text are shown by the dotted 
lines. 

Fig. 5 Momentum distribution of protons from the Lund string model for different (anti) neutrino 
beam energies (-Ebeam = 20 GeV - dashed line, -Ebeam = 50 GeV - solid line and -Ebeam 
= 80 GeV - dotted line. The experimental [^ cuts: W > 3GeV and x > 0.1 have been 
taken into account. The proton identification efficiency of the experiment is not included 
here. 

Fig. 6 Momentum distribution of protons from the quark-diquark model f^. The diquark 
fragmentation contribution is shown by the dashed line and the 4-quark fragmentation 
contribution by the dotted one. The details concerning the fragmentation function used 
are described in the text. 

Fig. 7 Momentum distributions of protons produced in the spectator mechanism in the neutrino 
(a) and antineutrino (b) induced deep inelastic scattering (solid line). The direct ttN 
contribution is shown separately by the dashed line, whereas the contribution from the 
valence quarks in the pion is shown by the dotted line. 

Fig. 8 Momentum distributions of protons corrected for the identification efficiency from the 
quark-diquark model (dotted line) and from the spectator model (dashed line) compared 
with the CERN experimental data [^. The solid line includes both mechanisms in the 
way described in the text. 

Fig. 9 The spectator mechanism for the production of slow pions in the (anti)neutrino induced 
reactions on the nucleon. 

Fig. 10 Momentum distribution of pions produced in the p(z/, 7r)X and p(v, 7r)X deep-inelastic 
processes. The predictions of the spectator model are shown as the solid [ttN) and the 
dashed (ttA) lines. The result of the Monte Carlo simulation according to the Lund string 
model P] are shown by the solid line with the centered points. 
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